Triplanar Model for the Gap and Penetration Depth in YBCO 
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YBaCuC>7 is a trilayer material with a unit cell consisting of a Cu02 bilayer with a CuO plane of 
chains in between. Starting with a model of isolated planes coupled through a transverse matrix 
element, we consider the possibility of intra as well as interplane pairing within a nearly antiferro- 
magnetic Fermi liquid model. Solutions of a set of three coupled BCS equations for the gap exhibit 
orthorhombic symmetry with s- as well as d-wave contributions. The temperature dependence and 
a-b in plane anisotropy of the resulting penetration depth is discussed and compared with experi- 
ment. 

preprint: cond-mat/9612106 
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I. INTRODUCTION 

Within the unit cell of YP^CuaCv, there are two 
Cu02 planes and a single set of CuO chains oriented 
along the 6-axis. The chains break the tetragonal sym- 
metry of the system and make it orthorhombic. The role 
of the chains is not yet fully understood but it is clear the 
they lead to a large in-plane anisotropy in various proper- 
ties such as the value of the zero temperature penetration 
depth and the DC conductivity etc. For example, in 
optimally doped YBa2Cu306.95, the penetration depth 
along the a-axis at T = is A„(0) ~ 1600A while that 
along the 6-axis is A b (0) - 1030A. @ At 300 K the ra- 
tio of the DC resistivities Pb/p a is larger than a factor of 
two. These large anisotropics would not exist in a purely 
tetragonal system. The role of the chains in YBCO has 
been investigated before usinga two plane model Q as 
well as a plane-chain model. ||[l0) While these calcula- 
tions have given us considerable insight, such a model is 
not fully realistic. Here we want to explicitly introduce 
the trilayer nature of the system. We start with a pair 
of isolated Cu02 planes and a set of chains (CuO plane). 
They are assumed to be coupled through a transverse ma- 
trix element, t\, between the chain and each of the planes 
and ti between the two planes. A set of three coupled BCS 
equations are used to describe the pairing with intra- as 
well as interband coupling. For the pairing potential, we 
use the form of the magnetic susceptibility introduced in 
the phenomenological nearly antiferromagnetic Fermi liq- 
uid model of Monien, Millis and Pines (mmp) which 
we multiply by a dimensionless parameter g a p which al- 
lows us to vary at will the strength of the diagonal (in 
plane) and off diagonal (interplane) couplings. Here the 
indices a[3 range over the three subbands. Numerical 
solutions of these equations using a fast Fourier trans- 
form technique give solutions which are an admixture of 
d, s and s symmetry. For the tetragonal crystal point 
group, these three symmetries would belong to two differ- 
ent irreducible representations and would, in general, not 
mix. But because of the existence of the chains, the sys- 
tem possesses reduced symmetry although, in most cases 



considered, nodes remain on the Fermi surface and the 
linear low temperature variation of the penetration depth 
remains. In section II, we give formalism. Solutions of 
the gap equations are discussed in section III along with 
results for the penetration depth. A brief conclusion is 
included in section IV. 



II. FORMALISM 

Consider a three plane per unit cell system with ener- 
gies £k,Q- Here a = 1, 2, 3 is an index that enumerates the 
three planes; two Cu02 which are identical to each other 
and one CuO plane consisting of chains oriented along the 
b-direction. The momentum index k is restricted to the 
two dimensional Cu02 Brillouin zone which has tetrag- 
onal symmetry. A tight binding model will be used to 
describe these isolated planes of the form: fll2|- 



2t a [(1 + e a ) cos(k x ) + (1 - e a ) cos(k y ) 
cos(fc x ) cos(fcy) - (2 - 2B a - fi a )] , 



- 2B, 



(1) 



where t a is the first neighbour hopping parameter, e Q an 
orthorhombic distortion which can be used to describe 
the chains, B a is the second neighbour hopping in units 
of t a , the ki are in units of reciprocal lattice vectors and 
Ha is related to the chemical potential and determines the 
filling. Here it will be treated as a parameter. In what 
follows, the taS will be taken to be the same for simplic- 
ity and denoted by t - The three dimensional nature of 
the system is introduced through transverse matrix ele- 
ments, t±, between the Cu02 planes and the CuO chains 
and £2 directly between the Cu02 bilayer. This changes 
the energies £k,a to band energies £k, a which we then 
use in the BCS gap equations assuming singlet pairing 
between equal and opposite momentum pairs. Including 
intraband (on diagonal) as well as interband (off diago- 
nal) pairing through the 2-body potential V^y , a p which 
describes scattering from k to k' with a and (3 band 
indices the gap in the a th band A^, a satisfies the BCS 
equations: pi 



1 
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^k',a 



k',/3 



2E ] 



tanh 



k'a 



\2kaTj' 



(2) 



In equation (g), f2 is the crystal volume, fee Boltz- 
mann's constant, T temperature and 



E u 



(3) 



are the quasiparticle energies of the a th band in the su- 
perconducting state. The single particle normal state 
band energies for a two layer model come from the Hamil- 
tonian: 



H = 



£i h 
t\ e 2 



with eigenvalues: 



£1,2 



£l + S 2 



± 



£l - £2 



+ *?, 



(4) 



(5) 



2 / V V 2 

and for the three layer model come from the Hamiltonian 



H 



£l tl t 2 
t\ e 2 tl 
t 2 ti E\ 



(6) 



with eigenvalues: 



£1 
£2 



£1 - t 2 

£1 + £2 + ti 



+ |\A £l ~ £ 2) 2 + 2t 2 {ei + £2 + t 2 ) + 8t 2 



and 



£3 = 



£1 + £2 + h 



e 2 ) 2 + 2t 2 (e 1 +e 2 +t 2 ) + 8t 2 



which are respectively the solid, dashed and dotted 
curves in Fig. ^. 

The set of three coupled equations (||) for the gap's 
Z\k, Q in the a th band (a — 1,2,3) are explicit once the 
pairing potentials 14:, k', a/3 are specified. Different models 
could be used. Here for simplicity and to be specific, 
we will assume that the superconducting condensation 
proceeds through the spin susceptibility Xk,k' which for 
the nearly antiferromagnctic Fermi liquid is peaked at 
the commensurate wave vector (71", tt) at the corner of the 
Brillouin zone. So as not to introduce new parameters, 
here we will use the form of Xk.k' introduced through a 
fit to NMR data by Millis, Monien and Pines (mmp). Q 
Their susceptibility and consequent pairing potential has 
the form: 



-t c 



k',a/3 = g a 



i + e 2 ik-k'-Qi 



(7) 



where the scale on He. k'. a/3 has been set by t for conve- 
nience and g a p are dimensionless parameters that deter- 
mine the overall size of the pairing as well as the relative 
amount of diagonal (a = [3) and off diagonal (a ^ (3) cou- 
pling desired. In all our calculations, the relative sizes of 
the g a j3 will be chosen at will with absolute value adjusted 
to get a critical temperature of 100K typical for the oxide 
superconductors. The other parameters in (0) are also 
specified and fixed at values given by MMP. £ = 2.5A is 
the magnetic coherence length and Q the commensurate 
wave vector (7r, tt) in the 2 dimensional CuC>2 Brillouin 
zone. 

In writing down the set of equations ([^), we have as- 
sumed that the coherence length transverse to the planes 
is short so that the gap can have different values in the 
different bands ek,a- The model still allows for a single 
gap, applicable to all three planes, to exist. The appro- 
priate limit is g a p independent of a and (3 in which case 
Z\k a on the left hand side of @ is clearly independent of 
a and a single gap parameter describes the entire system. 

Numerical solution of the set of equations (||) using fast 
Fourier transforms makes no assumption on the symme- 
try of the resulting gap parameters Z\k.a as a function of 
momentum k in the 2 dimensional Brillouin zone. In our 
work, we find that these solutions consist of a mixture of 
d x 2_ y 2 1 s and s x -2 +y 2 symmetries with higher harmonics 
as well as the more familiar lowest ones which are 



(s ) , 

s ° ~ Vk = 1 
s x 2 + y 2 ~ Vk = cos ( k x) + cos(k v ) 



"-x 2 —y 2 ~ ^?k 



l k = cos(k x ) — cos(fc !/ ). 



(8a) 
(8b) 
(8c) 



We are now ready to present results including the tem- 
perature dependence of the penetration depth A(T) for 
various model parameters. To obtain the penetration 
depths, we have adopted the following approximate pro- 
cedure. The lowest harmonics represented by (|8a| ) to ( [Sc] ) 
are projected out of our numerical solutions and these are 
used in all further calculations. This neglects higher har- 
monic contributions but is an adequate procedure since 
we are mainly interested in qualitative results and trends. 
In any case, our model for the pairing interaction (^) is 
simplified and is not believed to be accurate. We use a 
BCS approach while it is necessary to go to an Eliashberg 
formulation to get quantitative results. This goes beyond 
the scope of this work. 

For completeness, the approximate formula for the 
penetration depth employed here is: 



\-2 _ 



Aire 2 1 



c 2 fl^ 

k.a 



df(ik, a ) df(E Ka ) 



dEk 



(9) 



where the constant c is the speed of light, v a .i is the i th 
component of the electron velocity in the a th band 

1 <9e k , Q 



h dki 



(10) 
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and f(x) is the Fermi-Dirac distribution given by 

f(x) = (e^/^ + l)" 1 (11) 
We are now in a position to present results. 



III. RESULTS 

Before presenting results for the trilayer model, it is in- 
structive to consider a bilayer consisting of a Cu02 plane 
with tetragonal symmetry and a set of CuO chains di- 
rected along the b axis which has lower orthorhombic 
symmetry. This difference in symmetry has a drastic ef- 
fect on the solution of the BCS gap equation. To illustrate 
this, we take a model band structure: 



50.0 0.0 0.45 0.51 
50.0 1.0 0.00 1.2 



(12) 



with t in meV and the other parameters are dimension- 
less. We also consider the pairing potential to be of the 
mmp type and of equal magnitude in both chains and 
plane. So as to have a single critical temperature, we 
also include a smaller off diagonal pairing g a p ^ for 
Specifically, we take: 



9ap 



38.75 
1.0 



1.0 
38.75 



(13) 



Solutions of the BCS equations using fast Fourier trans- 
forms give an admixture of s a , s x 2 +y 2 and d x 2_ y 2 which 
all belong to the same irreducible representation for an 
orthorhombic system. The gap amplitudes are 



-0.196 
-0.644 



-0.467 9.482 
-0.182 1.482 



(d) 



(14) 



with the ^ as given in (| 

CuC>2 plane, the d x 2_ y 2 component is dominant with a 
small admixture of s and s x 2 +y 2 having the opposite 
sign (out of phase by it) . Although the same mmp inter- 
action as used in the CuC>2 planes has been assumed to 
also apply on the CuO chains, the very different symme- 
try of the quasi one-dimensional Fermi surface has led to 
a gap on the chains which is smaller by an order of mag- 
nitude. The d x 2_ y 2 component has the same phase as it 
has in the planes but it is now not much larger than its 
s and s x 2 +y 2 parts. In Fig. |](b) (middle upper frame) 
and (e) (middle lower frame), we show the Cu02 and 
CuO Fermi surfaces (dashed lines), respectively, as well 
as the contour of gap zeros (solid lines). In Fig. |](b), the 
Fermi surface has tetragonal symmetry but the gap ze- 
ros have nevertheless moved off the main diagonals (solid 
curve) because of the admixture into the mainly d x 2_ y 2 
gap function of s and s x 2 +y 2 components. As these are 



We note that in the 



negative in equation ([l4j), the zeros cross the Fermi sur- 
face at an angle (9) measured clockwise from the negative 
k y axis which is slightly larger than 45°. This is seen more 
clearly in Fig. 0(a) where we show the absolute value of 
the gap on the Fermi surface as a function of 9 for the 
plane. Note that the gap does not display full tetragonal 
symmetry in that the zeros are displaced from 45° and 
135° and that the gap at 90° is smaller than it is at 0° 
and 180°. 

For the chain case, the gap zeros (solid curves in Fig. 
|l|(e)) are displaced much more off the main diagonals 
and as illustrated in Fig. ||d) the value of the absolute 
value of the gap as a function of 9 along the quasi one- 
dimensional chain Fermi surface is very small. Note that 
this is a result of the orthorhombicity (ie, geometry of the 
Fermi surface) of the chain Fermi surface as the pairing 
potential has been assumed to have the same size and 
form in the chains as in the tetragonal planes. For the 
given value of pairing interaction, the gap ratio 2Z\ max /T c 
is 3.86 in the plane while it has been reduced to 0.47 in 
the chains. Finally, it is important to note that in Fig. 
|l|(d) , the gap on the Fermi surface starts at a finite value 
of 9 reflecting the geometry of the nearly linear chain 
Fermi surface. The very small energy scale exhibited by 
the gap in the chains leads directly to a very steep slope 
in the resulting 6-axis penetration depth as seen in the 
solid curve of Fig. |l|(f). The contribution of the chains 
to the a-axis penetration depth is by contrast (dashed 
line) very small. This is quite different from the contri- 
bution of the plane to the penetration depth illustrated 
in Fig. ^(c). In this case, the nearly tetragonal symmetry 
of the Fermi surface leads to very similar contributions 
for the temperature dependence of the a- (dashed) and b- 
(solid) direction penetration depth. The small differences 
between these two curves as to absolute value of A~ 2 (0) 
at zero temperature is due to the orthorhombic symme- 
try of the dispersion (equation (|l])). Unfortunately, these 
small differences are not measured directly in an experi- 
ment which probes the sum of contributions from chains 
and planes, ie, Fig. [j](f) and |l|(c). This sum is shown 
in Fig. ^ where it is compared with experiments (cir- 
cles). The agreement is poor although the order of the 
zero temperature slopes is correct, i.e. it is steeper in 
the chain direction (6-direction). The large discrepancy 
can be traced directly to the small energy scale found in 
the chains in this model and we can conclude that the 
data favours a gap which is of the same magnitude in 
the chains as it is in the planes. Finally, we note that 
A Q /A& ~ 1.4 in our calculation, which is reasonably close 
to (but not exactly) the measured value. 

It is easy to modify our model so that chains and planes 
have the same gap value. A single gap results for Z\k. a 
independent of a if all diagonal and off diagonal elements 
of the pairing potential are taken to have the same value, 
i.e. g a p = constant for any a and (3. In Fig. ^, we 
illustrate results obtained with the same band parameters 
as used in our previous model but now with g a p = 24.7 
for all a and (3 in which case a single gap results which 
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has the form: 



A k = [ -1.83 -0.288 10.09 ] 



9af3 



(15) 



This gap has mainly d-wave character with a small 
admixture of s and s x 2 +y 2 so that the gap zeros (solid 
curve) in Fig. |](b) and 0(e) have moved off the main diag- 
onal as in the previous example. Now these zero contours 
are the same in both chains and planes. The absolute 
value of the gap (on the Fermi surface) as a function of 
9 are shown in Fig. [|(a) and (d) for plane and chain, re- 
spectively, and the maximum value of 2A/k^,T c is 4.5 in 
the planes and 3.76 in the chains, much more comparable 
in value than was the case in the first model considered. 
It is to be noted that over the entire Brillouin zone the 
gaps in the planes and chains are not different but, of 
course, the Fermi surface geometry is totally different in 
these two cases. In particular, we note that on the Fermi 
surface the gap in the chains starts at a finite angle 9 near 
45° and that there is no Fermi surface at small values of 
9. 

The smaller value of the maximum gap found in the 
chains as opposed to the planes will lead to a steeper slope 
of the low temperature dependence of the inverse square 
of the penetration depth in the chain direction as seen in 
Fig. |^ where we compare with experiments. While theory 
predicts somewhat steeper slopes, qualitative agreement 
with experiment is obtained in this simple model. 

We turn next to results obtained for the trilayer model 
with CuC>2 planes on both sides of a CuO chain plane. 
The transverse coupling between CuC>2 and CuO plane 
is denoted by t\ and ti is a direct transverse coupling 
between the two Cu02 planes of the same unit cell. For 
{to,ti,t 2 } — {50, 30, 15} in meV the resulting Fermi sur- 
faces obtained from equation^ are shown in Fig. |5| where: 



{ to B a fl a } 




(16) 



with t in meV. The chain Fermi surface, denoted by the 
dashed line, is quasi one dimensional and has only or- 
thorhombic symmetry. The two Cu02 plane Fermi sur- 
faces are shown as dotted (odd) and solid (even) curves. 
Only the even band Fermi surface has tetragonal sym- 
metry while the odd band is flatter along k x than it is 
along k y . This orthorhombicity will lead to a different 
a- and ^-direction contributions to the penetration depth 
from the Cu02 planes themselves and this will now be 
described. 

A first set of results that is instructive is illustrated 
in Fig. |^. The chains are assumed not to contribute 
directly to the penetration depth but only provide an 
agency whereby the Cu02 plane acquires orthorhombic 
symmetry. This is accomplished easily if no pairing is as- 
sumed to act on the chains be it diagonal or off diagonal. 
That is, we take for the pairing potential the matrix: 



36.8 0.0 15.0 
0.0 0.0 0.0 
15.0 0.0 36.8 



(17) 



What is critical in equation ( |l7j ) is that the second row 
is all zero so the chains remain normal by arrangement. 
In this case, we obtain a gap of the form: 



1.81 -1.46 9.96 
0.0 0.0 0.0 
0.59 -0.80 10.88 



(«o) 

(d) 



(18) 



The gap on even (top row) and odd (bottom row) band 
is mainly d-wave like with s-components an order of mag- 
nitude smaller. The contours of gap zeros are illustrated 
as the solid curves in Fig. ||(b) and (e) for odd and even 
bands, respectively. The Fermi contour is the dashed 
curve in each case. The absolute value of the gap (\A\) 
as a function of angle 9 along the Fermi surface are given 
in Fig. [5](a) and (d) for odd and even band, respectively. 
In both cases, the gap zeros occur slightly before 9 = 45° 
and beyond 135° so that the gap does not have tetrag- 
onal symmetry even in the even band with a tetragonal 
Fermi surface. Also, the maximum gap on the Fermi 
surface is not quite the same in the two bands and is 
2Ajk^,T c =4.5 and 3.69 for even and odd cases, respec- 
tively. While these values are not very different, the con- 
tributions of even and odd bands to the a- and 6-direction 
penetration depths are very different. The odd band with 
orthorhombic Fermi surface gives a factor of 2 difference 
in the value of the zero temperature penetration depth as 
shown in Fig. ^(c). The solid curve applies along the b- 
direction and the dashed curve along the a-direction. For 
the even Fermi surface, which has tetragonal symmetry, 
the situation is very different as illustrated in Fig. ||(f) 
where it is seen that solid (along 6-direction) and dashed 
curve (along a-direction) are nearly the same (the slight 
difference is due to the difference in the o and b lattice 
parameters). Finally, in Fig. we compare results for 
the a- (dashed) and 6-direction (solid) penetration depth 
with experimental results. Again the curves are correctly 
ordered as to value of slope at zero temperature and the 
agreement is qualitative. It is as good as that shown in 
Fig. ||, although the physics of the two cases is quite dif- 
ferent. In Fig. 0, the chains are normal by assumption 
and do not contribute directly to the condensate density 
in the a- or 6-directions. The in plane anisotropy is due 
mainly to the odd band Fermi surface which has been 
made orthorhombic through its coupling to the chains 
and much less importantly to the s admixture in the 
gap function. Thus, we are dealing mainly with a band 
structure effect which dominates and overwhelms any gap 
symmetry effects. 

As a second illustration of a triplanar model, we show 
in Fig. ^ results for model parameters chosen so that 
the chains also become superconducting but by virtue of 
off diagonal coupling and not because there is intrinsic 
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pairing in the chains themselves. That is 522 is assumed 
to be zero. Specifically, we assume: 



5a/3 



44.4 10.0 0.0 
10.0 0.0 10.0 
0.0 10.0 44.4 



(19) 



while the band parameters remain the same as in the 
previous cases. In Fig. ||(b), (e) and (h), we show, re- 
spectively, odd, chain and even band Fermi surface as 
the dashed curves and contours of gap zeros is the solid 
curves. The gap solution is: 



A* 



1.21 -1.35 7.77 
0.65 -0.34 3.79 
-0.59 0.15 11.17 



M 
Vk 

(«») 

(d) 



(20) 



The absolute values of the gaps as a function of angle 
8 along odd, chain and even Fermi surface are shown in 
fig. 8(a), (d) and (g). Because no pairing was assumed 
to act in the chains, the gap is smallest in frame (d) with 
maximum 2Ajk^,T c equal to 1.88 as opposed to 4.58 and 
2.98 for even and odd band, respectively. Frames (c), 
(f) and (1) show that odd orthorhombic band makes a 
very different contribution to the 6-direction (solid curve) 
penetration depth than to the a-direction; the difference 
in zero temperature value being of order 2. The differ- 
ences are even more pronounced in the case of the chains 
which contribute very little to the a-direction (dashed 
curve) while the even tetragonal Fermi surface makes al- 
most equal contribution to a- and 6-directions. Because 
of the small energy scale present on the chains, which can 
be traced to the assumption <? 2 2 = 0, i.e. zero pairing on 
the chains, this model will not lead to good agreement 
with experiment. To achieve this, a final model is con- 
sidered with the same single gap for both CuC>2 planes 
and for the chains. If we take all g a p = constant = 18.4, 
we obtain a single gap 



Ak = [ -0.363 -0.502 10.14 ] 



Vk 

(»x) 

Vk 

(d) 

Vk 



(21) 



Detailed results are illustrated in the nine frames of 
Fig. [| In Fig. ^ (b), (e) and (h), the gap zeros are all 
the same (solid curves) yet because the Fermi contours 
(dashed curves) are different, the gap as a function of an- 
gle along the Fermi surfaces are different as shown in Fig. 
(a), (d) and (g) for odd, chain and even Fermi surfaces, 
respectively. The maxima in 2A/k^T c are respectively 
4.39, 4.55 and 4.07. Resulting contributions to the a- 
(dashed) and b- (solid) direction penetration depth are 
found in Fig. ^ (c), (f) and (i). The odd orthorhombic 
Fermi surface again contributes about twice as much to 
b- and to a- direction at low temperature while the even 
one contributes equally. The chains mainly contribute to 
the 6-direction. The sum of all contributions arc shown in 



satisfactory with a-direction curve (dashed) above the in- 
direction (solid) as observed. The ratio A /A& was ~ 1.6 
in the calculations, close to the observed value. 



IV. CONCLUSIONS 

We have presented results of calculations for the tem- 
perature dependence of the in plane penetration depth 
in a trilayer model of YBa2Cu307 which includes two 
CuC>2 planes and a set of CuO chains oriented along the 
6-direction. Transverse hopping matrix elements couple 
the three planes together leading to three coupled bands 
in which the gap parameters is an admixture of d x 2_ y 2 1 

S r 2 



x +y 



2 and s symmetries. While the s and d symmetry 
types belong to two different irreducible representations 
of the tetragonal crystal point group, they belong to the 
same representation in the orthorhombic case which is 
the case we are dealing with here because of the pres- 
ence of the chains. Nevertheless, it is found that the 
penetration depth remains linear in temperature at low 
temperature as a result of nodes in the gap crossing the 
Fermi contours. To compare with experiment, the ratio 
of A /A(, was kept near the observed value in all our cal- 
culations. With this constraint, it was found that the 
observed normalized temperature variation of the pene- 
tration depth in a- and 6-directions is difficult to under- 
stand in models where no intraplane pairing is introduced 
in the chains even if a large amount of interplane pairing 
is considered. The most natural explanation of the data 
is obtained when the pairing interaction is assumed not 
to depend significantly on sub-band and is taken to be 
independent of a and /?, i.e. Vk. k', up = Hc,k'- Also the 
zero temperature anisotropy in A a ,6 can be traced to the 
orthrhombicity of the band structure which resides in the 
chains and hybridized planes. 
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FIG. 4. Total contribution to the a- (dashed) and 
6-directions (solid) as a function of temperature to the in- 
verse square of the penetration depth normalized to its zero 
temperature value. The dots are experimental values Q with 
the 6-direction falling below the a-direction. The agreement 
is good and value of \ a /\b ~ 1.4 is also in reasonable but not 
exact agreement with experiment. 
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FIG. 5. The three Fermi surfaces in the triplanar model 
of YBCO. The chain like band (dashed curve) is nearly one 
dimensional, the odd band Fermi surface has orthorhombic 
symmetry (dotted curve) due to its interaction with the chains 
and the even band (solid curve) remains tetragonal. 
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FIG. 7. Total contributions to the inverse square of the 
penetration depth normalized to its zero temperature value 
as a function of temperature. The solid curve is for the 
6-direction and the dashed for the a-direction. The dots are 
the experimental data jl| with 6-direction curve below the 
a-direction curve. Results are for the trilayer model of Figs. 
|and§ in which it is assumed that there is no pairing on the 
chains on or off diagonal so they remain normal and do not 
contribute directly to the condensate. 
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FIG. 1. (a) and (d) are the absolute value of gap as a function of angle 8 measured clockwise from the k x axis along the 
Fermi surface for the Cu02 plane and CuO chains, respectively, (b) and (e) are Fermi surface contours (dashed line) and gap 
zeros (solid line) while (c) and (f) are the contributions to the inverse square of the penetration depth X~ 2 (T) as a function of 
temperature along a- (dashed curve) and 6-direction (solid curve). Note that the chains Fig. (d), (e) and (f) contribute almost 
exclusively to the 6-direction. The same pairing was included in chains and plane (equation H3) yet note the very small energy 
scale for the gap in the chains (d) as compare to the gap in the plane (a). This is due to the quasi-one dimensional nature 
of the chain Fermi surface. Note also that the low energy scale in the chains leads to the sharp drop in penetration depth in 
6-direction with increasing temperature seen in (f). 
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FIG. 3. Frame (a), (b) and (c) apply to the Cu02 planes while (d), (e) and (f) apply to the CuO chains with a nearly 
one dimensional Fermi surface. The pairing interaction is assumed to be the same on chains and planes and off diagonal so 
that a single gap results. Figs, (b) and (e) give the Fermi surface contours (dashed curve) and the gap zeros (solid curves) 
which are the same in both frames but different from pure d-wave which would be lines on the main diagonals in the Brillouin 
zone. Frame (a) and (d) give the absolute value of the gap as a function of angle 6 along the Fermi surface measured clockwise 
from the k x axis. Frame (c) and (f) give partial contributions to the inverse square of the penetration depth as a function of 
temperature along o-direction (dashed) and ^-direction (solid) . Note that the zeros in the gap are not at 45° and 135° as they 
would be in pure d-wave (Fig. (a) and (d)). 
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FIG. 6. Results for the triplanar model illustrated in Fig. q assuming no pairing diagonal or off diagonal on the chains 
which then remain normal. Frame (a), (b) and (c) apply to the odd band (orthorhombic symmetry) while (d), (e) and (f) apply 
to the even band (tetragonal symmetry). The first two frames (a) and (d) illustrate the absolute value of the gap along the 
Fermi surface as a function of angle 6 measured clockwise from the k x axis. These curves do not have full tetragonal symmetry 
and the zeros are displaced somewhat from the ideal d x 2_ y 2 case which would be at 45° and 135°. Frame (b) and (e) give the 
Fermi surface (dashed line) and gap nodes (solid line) for odd and even plane (tetragonal), respectively. The final frames (c) 
and (f) give partial contributions to the inverse square of the penetration depth as a function of temperature for a- (dashed) 
and b- (solid) direction. Note that the odd (orthorhombic) band gives a contribution to the 6-direction which is about twice 
its contribution to the a-direction, while the even band (tetragonal) gives nearly the same contribution in both directions. In 
the model, the ratio \ a /\b ~ 1.3. 
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FIG. 8. Results for a triplanar model with some off diagonal pairing on the chains but with 522 = 0, i.e. no direct pairing 
in the chains themselves. The chains become superconducting but with a small gap amplitude. Top frames apply to the odd 
(orthorhombic) Cu02 band, the middle three to the chain (CuO) band and the lower three frames to the even (tetragonal) 
Cu02 band. Frames (a), (d) and (g) give the absolute value of the gap as a function of angle 6 along the Fermi surface measured 
clockwise from the k x axis. Note the lack of tetragonal symmetry in the gap and the shift off 45° and 135° of the gap nodes. 
Frame (b), (e) and (h) give the Fermi surface (dashed curve) and gap zero contours (solid curve). The odd CUO2 plane band 
is orthorhombic as are the chains while the even CuCh band Fermi surface retains tetragonal symmetry. Contributions to the 
inverse square of the penetration depth as a function of temperature are given in frames (c), (f) and (i) for a-direction (dashed) 
and b-direction (solid). The odd orthorhombic band contributes almost twice as much to the 6- as to the a-direction while the 
even tetragonal band makes almost equal contributions. The chains contribute mainly in the b- direction. 
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FIG. 9. Results for a trilayer model with pairing interaction having the same value on and off diagonal, i.e. g a p — constant 
for all a and (3. This assumption ensures that a single gap applies to both the CuC>2 planes and the CuO chain. The three 
upper frames apply to the odd Cu02 orthorhombic band, the middle three to the CuO chain band and the lower three to the 
even Cu02 (tetragonal) band. Frames (a), (d) and (g) gives the gap on the Fermi surface as a function of angle 9 measured 
clockwise from the fc^-axis. Note the lack of tetragonal symmetry and the "missing parts" in the chain case which reflects the 
geometry of the quasi one dimensional Fermi surface. Frame (b), (e) and (h) give the Fermi contours (dashed line) for odd 
Cu02 (orthorhombic) band, chain like CuO band and even Cu02 (tetragonal) band. The solid curve are the gap zero contours 
and are the same in all three figures. Frames (c), (f) and (i) give the contributions to the inverse square of the penetration 
depth as a function of temperature for the a-direction (dashed) and the 6-direction (solid) . Note that the odd orthorhombic 
band contributes almost twice as much to the b- as it does to the a-direction in contrast to the even (tetragonal) band which 
makes nearly equal contributions. The chains contribute mainly to the b-direction. In this model part of the a-b asymmetry in 
the penetration depth comes from the orthorhombic nature of the even plane Fermi surface of frame (b). 
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FIG. 10. Total contributions to the inverse square ol the 
penetration depth normalized to its zero temperature value 
as a function of temperature. The solid curve applies to the 
6-direction and dashed to the a-direction. The dots are exper- 
imental data [jjj with the ^-direction falling below a-direction. 
The model is that of Fig. [] with pairing interaction assumed 
constant (i.e. g a p = constant) independent of a and /3 which 
leads to a single gap for all three planes in the primitive cell. 
This does not imply that it has d x 2_ y 2 symmetry because the 
bands do not have tetragonal symmetry and the gap contains 
an admixture of s and s x 2 +y 2 although it is still dominantly 

^x 2 — y 2 ■ 



